Self-organized criticality in single neuron excitability 
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We present experimental and theoretical arguments, at the single neuron level, suggesting that 
neuronal response fluctuations reflect a process that positions the neuron near a transition point 
that separates excitable and unexcitable phases. This view is supported by the dynamical properties 
of the system as observed in experiments, as well as by a theoretical mapping between the constructs 
of self organized criticality and membrane excitability biophysics. 



Cellular excitability is a fundamental physiological 
process whereby voltage-dependent changes in excit- 
ing and restoring membrane ionic conductances lead to 
an action potential (AP), a transient change in trans- 
membrane voltage. Hodgkin and Huxley [1 formalized 
a generic biophysical mechanism underlying the ignition 
and propagation of action potentials. In this formalism, 
as well as in its later extensions, the flow of ions down 
their electrochemical gradients is modulated by the prob- 
ability of ion channel proteins to reside in a conductive 
state. An extensive set of observations shows that the 
activity and response properties of neurons are highly 
variable, fluctuating over extended time scales in a com- 
plex manner (e.g. [2-4 .). Several approaches have been 
suggested to explain these fluctuations, largely focusing 
on the stochastic nature of underlying mechanisms [SHE] , 
non linearity and chaotic dynamics [TD] , or network 
level effects Here we examine variability and com- 
plexity in single neuron activity from a different angle, 
as reflecting self-organized critical fluctuations at a phase 
transition of excitability. Self-organized criticality (SOC) 
was suggested by Beggs and Plenz as a framework that 
explains various complex phenomena in neural systems 
at the network level [J2| [13] ; to the best of our knowl- 
edge application to single neuron excitability was never 
considered. 

Excitability is a lumped product of the individual 
states of numerous interacting ion channels. Focusing 
on changes in excitability over time far beyond the mil- 
lisecond scale of a single action potential, we define the 
excitability status of a cell as the aggregated, macroscopic 
availability of ionic channels to move into the conduc- 
tive state and participate in the generation of action po- 
tentials. Defined as such, the excitability status of a 
neuron reflects the susceptibility of the cell to produce 
an action potential in response to input above a given 
amplitude. It is instructive to think of excitability sta- 
tus in the context of G max , the maximal conductance 
(or number of ionic channels) in a unit area of mem- 
brane; in the original Hodgkin and Huxley formalism, 
aims at the scale of milliseconds, maximal conductance 
is a structural constant that sets limits on the instanta- 



neous input-output relations of the membrane. But when 
long-term effects arc sought, the maximal conductance 
might (and indeed should) be treated as a (macroscopic) 
system variable governed by the stochastic dynamics of 
ion channels, reflecting their activity-dependent tran- 
sitions between available and unavailable microscopic 
states. These microscopic transitions are globally and 
locally coupled via membrane voltage, ionic concentra- 
tions and cellular modulatory and homeostatic processes; 
their impacts on excitability status depend on the specific 
type of ionic channel involved (i.e., mediating exciting 
or restoring ionic flows). Excitability status translates 
into actual stimulus evoked response in a highly non- 




t 500 1000 1500 2000 2500 



lime (s) ^ms 

FIG. 1. Dynamics of neuronal excitability status. (A) The 
effect of modulating Gjv a in HH model. As excitability status 
decreases, the AP is delayed. Below a certain threshold, no 
AP is produced. (B) AP latency in (A) as a function of a, 
a factor modulating the nominal GNa,max in the HH model, 
demonstrating the existence of a sharp threshold. (C) Exper- 
imental demonstration of excitability dynamics. Following a 
sequential stimulation (ordered top to bottom), the AP is de- 
layed, reflecting a decrease in the excitability status of the 
neuron. (D) The AP latency plotted as a function of time in 
an experiment where the stimulation rate is changed. For low 
stimulation rate, the excitability status stabilizes at a fixed, 
supra threshold value. For high stimulation rate (20Hz), ex- 
citability status decreases below threshold, and the neuron 
responds intermittently. (E) The voltage traces of the 20Hz 
stimulation in (D), together with the following 1Hz stimula- 
tions. 
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linear threshold-governed manner: as excitability status 
decreases, the time delay from the stimulus to action po- 
tential is extended; below a threshold level of excitability 
status, no action potentials can be generated (FigurellK). 

Modelling the dynamics of excitability over extended 
time scales is not a trivial matter. On the one hand, the 
multitude of underlying processes renders it unlikely that 
low dimensional deterministic models account for tempo- 
ral long-range correlated fluctuations in excitability sta- 
tus [2]. On the other hand, stochastic models - while be- 
ing more natural given the known microscopic machinery 
- are limited in their capacity to account for the practi- 
cally unbounded temporal complexity and structure seen 
in the data [T3] . 

The concept of Self-Organized Criticality [H] desig- 
nates a cluster of physical phenomena characterizing sys- 
tems that reside near a phase transition. What makes 
SOC unique is the fact that residing near a phase tran- 
sition is not the result of a fine-tuned control parameter; 
rather, in SOC the system positions itself near a phase 
transition as a natural consequence of the underlying in- 
ternal dynamic process that pushes towards the critical 
value. Such systems exhibit many complex statistical 
and dynamical features that characterize behavior near a 
phase transition, without these features being sensitive to 
system parameters. Dickman and his colleagues [T71 [TS] 
formalized a scheme for generating SOC from a conven- 
tional system exhibiting a phase transition. They have 
shown that many of the canonical models of SOC, includ- 
ing sandpile and forest fire models, are in fact absorbing 
state (AS) systems, amended with a carefully designed 
feedback: dissipating energy whenever the system is su- 
percritical (i.e. permanently active without settling into 
AS), and driving the system whenever it is subcritical 
(i.e. when and only when it settles into the AS). 

Several theoretical and phenomenological observations 
bring us to consider the framework of SOC as a useful 
framework to account for excitability fluctuations. The 
macroscopic state of the membrane can be divided into 
two distinct phases, determining its response to a given 
input: an excitable phase and an unexcitable phase, sepa- 
rated by a sharp boundary in parameter space. These pa- 
rameters, however, are dynamic on slower scales, and are 
a function of the system state. The dynamics are driven 
by neural activity, which serves as a temperature-like pa- 
rameter, and the single AP serving as a drive (quantal 
influx of energy, or small increase in temperature) . How- 
ever, neural activity level is not a control parameter set 
by the experimenter; rather, it depends, in turn, on the 
excitability of the neuron, giving rise to the feedback loop 
inherent to SOC. In the absence of activity, the neuron 
reaches an excitable state, while increased activity re- 
duces the excitability status, and (when high enough) 
pushes the membrane into the unexcitable phase. Of 
course, not all classes of neurons follow this simplistic 
process, but the general idea holds - activity pushes ex- 
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FIG. 2. Steady state characterization of the response. (A) 
response latencies (solid line) in response to a stimulation se- 
quence with slowly increasing stimulation rate (dashed line). 
(B) Fail (no spike) probability as a function of stimulation 
rate. A critical stimulation rate is clearly evident. (C) Mean 
response latency as a function of stimulation rate. The in- 
crease of the latency accelerates as the stimulation rate ap- 
proaches the critical point. (D) The jitter (coefficient of vari- 
ation) of the latency as a function of stimulation rate. 



citability status towards a threshold state, while at the 
longer time scale regulatory feedback pulls the system 
back. 

In what follows, we briefly summarize experimental ob- 
servations on excitability dynamics in biological neurons 
exhibiting footprints of SOC, and discuss the mapping 
of membrane biophysics onto this framework. Fluctu- 
ations of excitability are monitored using extracellular 
recording and stimulation in cultured single neurons, iso- 
lated from their networks by means of pharmacological 
synaptic blockage [3]. The neurons are stimulated with 
sequences of short, identical electrical pulses. For each 
pulse, the binary response (action potential produced or 
not) is registered, as well as the latency from stimulation 
time to the voltage peak of the recorded action potential. 
The amplitude of the stimulating pulses is set such that 
the neuron will respond in a 1:1 manner under low rate 
(1Hz) conditions. When stimulation rate is abruptly in- 
creased to a higher value, the latency to the action poten- 
tial increases and stabilizes on a new value (Figure []J)). 
At a certain stimulation rate ro, the 1:1 response mode 
(the stable regime) breaks down, and the neuron starts 
to respond intermittently, with irregular spiking and jit- 
tered latency (the intermittent regime) . The steady state 
properties of the two response regimes may be observed 
by slowly changing the stimulation rate. As seen in the 
result of this 'adiabatic' experiment of Figure [2j the sta- 
ble regime is characterized by 1:1 response (no failures), 
stable latency (low jitter) and monotonous dependency 
of latency on stimulation rate. The intermittent regime 
is characterized by a failure rate which increases with 
stimulation rate, unstable latency (high jitter) and inde- 
pendence of the mean latency on the stimulation rate. 
Figure [2p exemplifies the role of the stimulation rate in 
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such experiments: it externally sets an upper limit to the 
activity rate of the neuron. The activity rate itself is a 
dynamic variable of the system. Figure [2] also indicates 
that change in both mean latency and its variance accel- 
erates as the critical stimulation rate is approached. The 
exact value of ro can vary considerably between neurons, 
but its existence is observed in practically all measured 
neurons (see details in [4]). 

When repeated stimuli at a fixed rate above ro are ap- 
plied, failures to induce action potentials undergo fluctu- 
ations characterized by scale-free long-memory statistics. 
The power spectral density (PSD) exhibits a power-law 
tail at the low frequency domain. The characteristic ex- 
ponents do not depend on the stimulation rate, as long 
as the latter is kept above ro (Figure [3]A_) . The typical 
exponent of the rate PSD is j3 = 1.26 ± 0.21 (mean ± 
SD, calculated over 16 neurons). Clearly, the approach 
to ro from below resembles an approach to a phase transi- 
tion; the critical characteristics, however, are maintained 
throughout the range above ro. 

Within the intermittent regime, the distributions of 
the lengths of consecutive response sequences (i.e. peri- 
ods of time the neuron is fully excitable, responding in 
a 1:1 mode) and consecutive no-response sequences (i.e. 
periods of time the neuron is not responding at all), are 
qualitatively different (Figures [3^3 and[3p). The consec- 
utive response sequence length histogram is strictly ex- 
ponential, having a characteristic length; the consecutive 
no-response sequence length histogram is wide, to the 
point of scale-freeness (power law distribution) , suggest- 
ing that the fluctuations are dominated by widely dis- 
tributed excursions into an unexcitable state. Moreover, 
as shown in Figure [3)3, during the intermittent regime 
the response of the neuron is characterized by switches 
between quasi-stable modes, typical temporal patterns 
that dominate the response sequence. 

Taken together, the experimental results summarized 
above support the interpretation of excitability fluctua- 
tions as reflecting a SOC state of the membrane. This in- 
terpretation may also be theoretically supported, at least 
to some extent, by considering the underlying biophysical 
machinery. The state of the membrane is a function of 
the individual states of a large population of interacting 
ion channel proteins. A single ion channel can undergo 
transformations between uniquely defined conformations, 
conventionally modeled as states in a Markov chain. The 
faster transition dynamics between states is the founda- 
tion of the Hodgkin Huxley model, which describes the 
excitation event itself - the action potential. But for the 
purpose of modeling the dynamics of excitability, rather 
than the generative dynamics of the action potential it- 
self, it is useful to group these conformations into two 
sets [lOl HH |20j : the available, in which channels can 
participate in generation of action potentials, and the 
unavailable, in which channels are deeply inactivated and 
are "out of the game" of action potential generation. The 




FIG. 3. Scale free fluctuations in the intermittent mode. (A) 
Periodograms of the failure rate fluctuations, at 5 different 
stimulation rates above ro. (B) Length distribution of spike- 
response sequences, on a semi logarithmic plot, demonstrating 
an exponential behaviour. Example from one neuron stimu- 
lated at 20Hz for 24 hours. (C) Length distribution of no-spike 
response sequences from the same neuron, on a double log- 
arithmic plot, demonstrating a power-law-like behavior. (D) 
Pattern modes in binary response sequences. Extracts (ap- 
proximately lOmin long) from the response pattern of two 
neuron to long 20Hz stimulation. White pixel represents a 
spike response, black represents no-spike response. The re- 
sponse sequence is wrapped. Transitions between pattern 
model are evident. 

microscopic details of the single channel dynamics in this 
state space, not to mention the collective dynamics of the 
interacting ensemble, are complex [2"lTf2"3"] and no satis- 
factory model exists to date. 

However, it has been suggested recently [101 [20] that 
the transition dynamics between the available and un- 
available states may be expressed in terms of an "adap- 
tive rate", Logistic-like model of the general form: 

x = -f{ 1 )x + g{x){l-x), (1) 

where / is a function of the neural activity measure 7, 
and g{x) is a monotonically increasing function of the 
system state x. 

Following the lead of the above adaptive rate approach, 
one can consider, for instance, a model in which x repre- 
sents the availability of a restoring (e.g. potassium) con- 
ductance, and the state of the single channel is repre- 
sented by a binary variable af, a-i — is the unavailable 
state and <7j = 1 is the available state. Unavailable chan- 
nels are recruited with a rate of x, while available chan- 
nels are lost with a rate of 2 — 7. This picture gives rise 
to a dynamical mean field like equation: 

x = (7 - l)x - x 2 . (2) 

This toy model is a variant of a globally coupled Contact 
Process (CP), a well-studied model exhibiting an absorb- 
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ing state phase transition [21]. Here, x = is the absorb- 
ing state, representing the excitable state of the system. 
In the artificial case of 7 as an externally modified control 
parameter, for 7 < 1 (low activity) the system will always 
settle into this state, and the neuron will sustain this level 
of activity. For 7 > 1, the system will settle on x* = 7— 1, 
an unexcitable state, and the neuron will not be able to 
sustain activity. Feedback is introduced into the system 
by specifying the state dependency of 7: An AP is fired 
if and only if the system is excitable (x=0), giving rise 
to a small increase in 7. When x > 0, the system is un- 
excitable, APs are not fired, and 7 is slowly decreased. 
This is an exact implementation of the scheme proposed 
in pTJ [18] : an AS system, where the control parameter 
(activity, 7) is modified by a feedback from the order pa- 
rameter (excitability, a function of x). As always with 
SOC, the distinction between order and control parame- 
ters becomes clear only when the conservative, open-loop 
version of the model is considered. Note that the natural 
dependency of the driving event (the AP) on the system 
state resolves a subtlety involved in SOC dynamics: the 
system must be driven slowly enough to allow the ab- 
sorbing state to be reached, before a new quantum of 
energy is invested. In most models, this condition is met 
by taking driving rate to be infinitesimally small. 

In summary, we have given several arguments, exper- 
imental and theoretical, in support of a connection be- 
tween the framework of SOC and the dynamics under- 
lying response fluctuations in single neurons. Naturally, 
the simple model leading to equation [2j is not unique, 
probably wrong in its microscopical details. Moreover, 
excitability is determined by more than one order param- 
eter, and the interaction types are much more heteroge- 
neous, controlled by an aggregate of such equations, rep- 
resenting the exciting and restoring forces, each pushing- 
pulling excitability to opposite directions. Nevertheless, 
while respecting the gap between theoretical models and 
biological reality, SOC seems to capture the core phe- 
nomenology of fluctuating neuronal excitability, and has 
a potential to enhance our understanding of physiological 
aspects of excitability dynamics. 
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